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Abstract. We construct a commutative algebra J^x of difference operators 
in W , depending on p + 3 parameters which is diagonalized by the multivari- 
able Racah polynomials Rp{n\x) considered by Tratnik |27) . It is shown that 
for specific values of the variables x = {x\,X2, ■ ■ ■ ,Xp) there is a hidden du- 
ality between n and x. Analytic continuation allows us to construct another 
commutative algebra An in the variables n = (ni, n2, • • • , Hp) which is also di- 
agonalized by Rp{n; x). Thus Rp{n; x) solve a multivariable discrete bispectral 
problem in the sense of Duistermaat and Griinbaum [8]. Since a change of the 
variables and the parameters in the Racah polynomials gives the multivari- 
able Wilson polynomials 1261 . this change of variables and parameters in Ax 
and An leads to bispectral commutative algebras for the multivariable Wilson 
polynomials. 
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1. Introduction 

It is well known that any family of polynomials orthogonal with respect to a 
positive measure supported on the real line satisfies a three term recurrence relation, 

di{n)pn+i{x) + do{n)pn{x) + d_i(n)p„_i(x) = xpn{x). (1.1) 

Furthermore it was shown by Bochner 5J that the classical orthogonal polyno- 
mials of Jacobi, Hermite and Laguerre can be characterized by the fact that they 
satisfy a second order differential equation of the form 

£ d 
C2{x)-^Pn{x) + ci{x)—pn{x) + co(a;)p„(x) = n{n)pn{x). (1.2) 

Thus the classical orthogonal polynomials are simultaneously eigenfunctions of a 
difference operator in n and a differential operator in x, i.e. they are part of a 
discrete-continuous version of the bispectral problem. The bispectral problem is 
concerned with finding differential or difference operators that have common eigen- 
functions. More precisely, find a differential or a difference operator L„, indepen- 
dent of X, acting on functions of n, and a differential or a difference operator 
which is independent of n acting on functions of x, so that 

Lnilj{x, n) = f{x)'il){n, x) 

Lxijj{x, n) = n{n)'ip{n, x), 

where f{x) and iJ-{n) arc functions of x and n respectively. This problem has 
appeared in many areas of mathematics, physics and engineering such as limited 
angle tomography, soliton equations and their master symmetries, particle systems, 
algebraic geometry and representations of infinite dimensional Lie algebras, see for 
instance [H [lOl [U [H H [T2] as well as the papers in p3]. 

While in (jl.ip and (|1.2|) L„ is a difference operator and Lx is a differential op- 
erator the bispectral problem considered in Duistermaat and Griinbaum [8] was 
for both Ln and to be differential operators which give rise to the continuous- 
continuous part of the bispectral problem. An advantage of this is that it puts 
both variables on an equal footing so it may be possible to obtain L„ from or 
vice versa via some hidden symmetry. In order to pursue this strategy for orthog- 
onal polynomials we look for more general families of polynomials which have the 
classical polynomials as limiting cases and for which L„ and are difference oper- 
ators. Indeed, for the discrete classical orthogonal polynomials (Charlier, Meixner, 
Krawtchouk and Hahn) equation (|1.2p is replaced by a difference equation of the 
form 

Ci(x)p„{x + 1) + Co{x)p„{x) + C^i{x)pnix - 1) = /i(n)p„(a;), (1.4) 
see for instance |23) . while (jl.ip is the same. The symmetry between n and x, which 
allows us to connect the two equations, is transparent in the cases of the Charlier, 
Meixner and Krawtchouk polynomials from their explicit formulas in terms of hy- 
pergeometric series. However, this is no longer true for the Hahn polynomials. One 
essential difference between equations (|1.1|) and (|1.4p for Hahn polynomials is that 
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the eigenvalue is quadratic in n, while the eigenvalue in (jl.ip is linear. Also 
the coefficients of the difference operator ()1.4p are polynomials of x, while the coef- 
ficients in (jl.ip are rational functions of n. Thus, we naturally arrive at the Racah 
polynomials, where the eigenvalue on the right-hand side of (|l.ip is also quadratic, 
and the symmetry between n and x can be seen from the 4^F^ representation. A 
similar duality holds also for the most general family - the Askey- Wilson polyno- 
mials [2 . In fact, the duality is even deeper and it can be extended beyond the 
polynomial case, to bi-infinite difference operators 13J. 

One beautiful extension of the above theory to orthogonal polynomials of more 
than one variable is related to the theory of symmetric functions known as Macdonald- 
Koornwindcr polynomials, see for instance [HI [201 EH Ell ■ A deeper understand- 
ing of the duality stems from affinc Hecke algebras considerations, introduced by 
Cherednik 7] in his proof of some conjectures about Macdonald polynomials (see 
also [23] for non-reduced root systems i?C„). 

In a different vein, one can also look for multivariable orthogonal polynomials 
satisfying an equation similar to (|1.4p . The first obstacle is the fact that in the case 
of more that one variable we have more than one monomial of total degree n, for 
n > 0. A natural classification of the possible weights should be independent of the 
way we order the monomials of a given degree when we apply the Gram-Schmidt 
process. This led to the notion of admissible difference operators, which was used 
by Iliev and Xu ^16] to characterize discrete multivariable orthogonal polynomials, 
satisfying (partial) second-order difference equations analogous to (ll.4|l . It is inter- 
esting to notice that the multivariable versions of Hahn, Krawtchouk and Meixner 
polynomials considered there had already appeared in the literature in different 
applications , 18, _22j . An intriguing question left open in [16| is whether these poly- 
nomials are bispectral and if the bispectrality can be derived from duality. Notice 
that this time a difference equation similar to (jl.ip is no longer an automatic con- 
sequence of the orthogonality, and its existence will depend on the way we order 
the monomials of total degree n. By analogy with the one dimensional theory if we 
hope to get the bispectrality from duality we need to go at least to multivariable 
Racah polynomials. 

In the present paper we answer afhrmatively the questions raised above by con- 
sidering the multivariable Racah polynomials defined by Tratnik '27] . These poly- 
nomials are orthogonal in W with respect to a weight depending on p -I- 2 real 
parameters /3o, • ■ ■ , Pp+i and a positive integer N . First, we define a difference 
operator Cp{x; f3; N) in the variables xi, X2, ■ ■ ■ , Xp which is triangular (see Defi- 
nition [2T3|) and self-adjoint with respect to the multivariable Racah weight. Here 
we have put /3 = (/3o, . . . , Pp+i) E Rp+^ and x = {xi,X2, ■ ■ ■ , Xp) G W. As an 
immediate corollary we see that any set of orthogonal polynomials with respect to 
Racah inner product will be eigenfunctions of Cp{x; /3; N). 

It is a bit surprising that this difference operator has a much more complicated 
structure than the operator for Hahn polynomials. Besides the fact that its coef- 
ficients are rational functions of x, it is supported on the cube { — 1,0,1}''. This 
means that it has 3^ rational coefficients, compared to p{p + 1) polynomial coeffi- 
cients for the operator corresponding to multivariable Hahn polynomials. In order 
to define and study the basic properties of this operator, we introduce p involutions 
on the associative algebra of difference operators in x with rational coefficients, see 
Section [2| for more details. Similar involutions (on the dual n-side) played a crucial 
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role in [111 I13j in the construction of orthogonal polynomials satisfying higher-order 
differential or g-difference equations, see Remark 14.81 for more details. 

Next we focus on the basis of orthogonal polynomials {Rp{Ti;x; p; N) : n G 
Nq such that |ri| = rii + • • ■ -|-np < N} constructed by Tratnik. We define a commu- 
tative algebra Ax generated by p difference operators in the variables xi,X2, ■ ■ ■ ,Xp, 
which is diagonalized by {_Rp(n; x; /3; A*")}. An analytic continuation argument 
shows that these equations are valid even when iV ^ Nq, in which case the polyno- 
mials {Rp{n] x; P; N)} are defined for all n e Nq. All this is the content of Section 

m 

In the next section, we prove that appropriately normalized, the polynomials 
{Rp{n; x; f3; N)} possess a certain duality between the variables n and x. The proof 
is obtained by applying twice an identity of Whipple connecting the terminating 
Saalschiitzian ^F^. The duality allows us to define an isomorphic commutative 
algebra An of difference operators in the variables ni,n2, ■ ■ ■ ,np which is also di- 
agonalized by {Rp{n; x\ /3; A'^)}, thus proving the bispectrality. 

In Section [5] we illustrate the bispectrality of other families of multivariable or- 
thogonal polynomials. In Subsection 15.11 we show that a change of the variables 
and the parameters in the Racah polynomials gives the multivariable Wilson poly- 
nomials ^26J . Therefore, changing the variables and the parameters in Ax and An 
leads to bispectral commutative algebras for the multivariable Wilson polynomials. 
In Subsections 15 . 2tf5 .31 we consider limiting cases which lead to bispectral commuta- 
tive algebras for multivariable Hahn and Jacobi polynomials. The explicit formulas 
derived there show that although the operators on the x side simplify tremendously 
when we go to the Hahn or Jacobi polynomials, the operators on the n side are more 
or less the same. In the case of Jacobi polynomials. Ax becomes a commutative 
algebra of differential operators going back to the work of Appell and Kampe de 
Feriet [1]. In the last subsection we illustrate the bispectrality and the duality of 
the multivariable Krawtchouk and Meixner polynomials. 

We have collected explicit formulas in dimension two in an Appendix. The reader 
eager to get a feeling about the complexity of the formulas or willing to look at 
some examples first is referred to the Appendix. 

2. Triangular difference operator in W 

In this section we introduce a triangular difference operator and we discuss its 
basic properties. 

2.1. Basic notations and definitions. Consider the ring V — M[a;i,X2, . . . ,Xp] 
of polynomials in p independent variables x\^xi^ . ■ ■ iXp with real coefficients. We 
denote by V\ = IR.[Ai, A2, . . . , Ap] the subring of V consisting of polynomials in the 
variables Ai, A2, . . . , Ap, where 

\i ^ \i{xi) = Xi{x.i ^ Pi) for I = 1, 2, . . . ,p, (2.1) 

and Pi, P2, ■ ■ ■ , Pp are parameters. 

For i = 1, 2, . . . ,p we define an automorphism li of V by 

= —Xi — Pi and Ii{xj) — Xj for j 7^ i. (2-2) 

Notice that li is an involution (i.e. holi = Id). The importance of these involutions 
is summarized in the following remark. 
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Remark 2.1. The involutions li preserve V\. Conversely, if a polynomial q ^ V 
is preserved by the involutions li for i — 1,2, . . . ,p then q G Vx- 

Let {ei, 62, ... , Bp} be the standard basis for W. We denote by Ex^, A^;. and 
Vxi, respectively, the customary shift, forward and backward difference operators 
acting on functions f of x — {xi,X2, ■ ■ ■ , Xp) as follows 

ExJix) = f{x + e,) 

A../(x) = fix + eO - fix) = iEx, - l)fix) 

V. = fix) - fix - e.) = (1 - 

Throughout the paper we use the standard multi-index notation. For instance, 
if = (i^i, 1^2, ■ ■ ■ , i^p) G ^ then 

X — x-^ X2 Xp , ^ -'^xi^x-i 

and = 1^1 + ^2 + • • • + I'p. 

Let I'j: be the associative algebra of difference operators of the form 

where S" is a finite subset of Z'' and l^ix) are rational functions of x. Thus, the alge- 
bra Vx is generated by rational functions of x, the shift operators Ex^ , Ex^ , ■ ■ ■ i Ex^ 
and their inverses E~^, E~^, . . . , E^^ subject to the relations 

Ex, ■ gix) = gix + ej)Exj, (2.3) 

for all rational functions gix) and for j — 1,2, . . . ,p. For every i £ {1,2, . . . ,p} the 
involution li can be naturally extended to Vx, by defining 

I'M^)) = gihix)), hiEx,) = E-\ hiEx,) = Ex, for j ^ i. (2.4) 

It is easy to see that li is correctly defined because the relations \2.i\ are preserved 
under the action of li, i.e. we have 

HEx,) ■ higix)) = higix + e,))hiEx,), 

for i,j G {1,2, . . . ,p}. 

For the sake of brevity we say that L e Vx is /-invariant if Ij (L) = L for all 
3 G {l,2,...,p}. 

2.2. The operator £p. Below we define an operator which is /-invariant. Clearly, 
an operator which is /-invariant is uniquely determined by its coefhcient of E'^ (or 
equivalently AJ!) with Vi > for i = 1,2, . . . ,p. 

Let 1/ = (i^i, 1/2, . . . , i^p) G {0, 1}'' \ {0}^, and let {vi^, Ui^,. . . , Vi^} be the nonzero 
components of v with < ii < ^2 < ■ • ■ < is < P + 1 and s > 1. Denote 
Au ^ixi^ -l-Ai -/3o)(a;ji +Ai) 

^ nfc=2('^'fc +^ifc-l +/3ifc)(^ifc + + + 1) ^2 

X (x,^ -f/3p+i + iV)(iV-a;iJ. 

An arbitrary v E 1,'^ can be decomposed as v = — , where G Ng with 

components z^^*" = max(zyj,0) and j/^^ = — min(jA,, 0). For v G { — 1, 0, 1}'° \ {0, 1}^ 
we define 

A,^ iA,+ +,-). (2.5b) 
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Here is the composition of the involutions corresponding to the positive coor- 
dinates of . 

Finally, we define the operator 



/:p-/:p(xi,...,a;p;/3o, /?!,..., /3p+i;iV)= V (-1)I''"U,A^ 



!^e{-i.o,i}p\{o}J' 

(2.6) 

In the above formula we use again multi-index notation, i.e. denotes the product 
Y^j^i^x] and similarly stands for the product O^'^i ■ Since /i(Aa;J — 
—Vxi and Ii{Axj) — for i ^ j, the operator defined by equation l|2.6p is 
/-invariant. 

Lemma 2.2. Let L G be an I -invariant difference operator. If 

p 

Y[{2x^ + l3^)L{q) e V for every q G Vx, 

i=l 

then L preserves V\, i.e. L{V\) C Vx. In particular, the operator Cp defined by 
T5l) - (|2.6p preserves Vx. 



Proof. Let q E Vx- Since both L and q are /-invariant, it follows that L(q) is also 
/-invariant. Therefore, it is enough to prove that L(q) G 'P, since Remark 12.11 then 
implies that L{q) G Vx- If we write L{q) as a ratio of two polynomials 

where qi and (72 have no common factors, then apart from a nonzero multiplicative 
constant, q2{x) must be the product Ilie-R'^^s;! + where K C {l,2,...,p}. 
Assume that /C 7^ and let i G K. Since Ii{2xi + f3i) = —{2xi + Pi), we see that 
IiiQ2ix)) = —q2{x) and therefore Ii(qi{x)) — —qi{x). But it is easy to see that 
this condition forces qi{x) to be an odd polynomial of {2xi + Pi), which shows that 
{2xi + Pi) is a common factor of qi and q2 leading to a contradiction. 

It remains to show that Cp satisfies the conditions of the Lemma. This follows 
easily from formulas (j2.5p - (|2.6p combined with the fact that Aa;. (Af ) is divisible by 
A^.(A,) = (2a;, + /3, + 1), and V^^iX^) is divisible by V^^iXi) = (2a;, -f - 1) for 
every fc G N. □ 

For d G No let us denote by "P^ the space of polynomials in Vx of (total) degree 
at most d in the variables Ai, A2, . . . , Ap, with the convention that V^^ = {0}. 

Definition 2.3. We say that a linear operator L on Vx is triangular if for every 
d G No there is c^; G M such that 

L{q) = Cdq mod V^'^ for all q G V^. 

The main result in this section is the following proposition. 

Proposition 2.4. The operator Cp defined by (|2.6p is triangular. More precisely, 
we have 

Cp{q) = -d{d - 1 + Pp+i - Po)q mod Vt^ for every qeVi (2.7) 
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Proof. For j £ {1, 2, . . . ,p} let 

i^e{-i,o,i}''\{o}'' 
where are given by (|2.5p . Then 

It is easy to see that for every j, the operator is /-invariant, satisfying the 
conditions of Lemma [2.2[ and therefore Cpj preserves P\. Moreover, since A^, is 
a ratio of a polynomial of degree 2(|i/"'"| + + 1) over a polynomial of degree 
2(1^^1 + and since Ax^ and Vx^ decrease the total degree of a polynomial by 

1, one can deduce that for j > 3 we have Cpj{Vf) C P'^^^. Thus 

Cp{q) = Cp.iiq) + Cp.2{q) mod P^^-^ for every q £ P^ (2.8) 
i.e. we can ignore the operators Cpj for j > 3. Clearly, it is enough to prove (|2.7[) 
for g = A" = A^i'^A^' • • • Ap^ Let us write Cp^i = YlLi 4!i' ^^^^^^ 

and 

, _ {X^ +|3^- Po){x, + P,){xi + + jV)(jV - Xj) 

(2x,+A)(2x.+A + l) ^^-^^^ 

. _ X^{X^ + /3o)(A^ - - A + /3p+l)(iV + + A) 

(2x, + /?,)(2x.+/3.-l) • ^ ^ 

Using the binomial formula one can check that 

/rii-l \ 



\ fc=0 / 

(2.10a) 

V.,(Ar') = (V..A,) ( £ Af(x,)Ar-^-'=(x, - 1)] 

\ k=0 ) 

(2.10b) 

From (12. 9p . (I2.10|) and the fact that the operator c!"^\ satisfies the conditions in 
Lemma 12.21 one can deduce that 

4a(^r') = K(/3o - + 1) - n?)xf ' + O(x^-i), 
which shows that 

(A") = A"— ^-^W (AD = (^.(/?o - /3p+i + 1) - ri^A" mod T^I-I^V 

Thus 

^P,i(A") = f]4:l(A") = ( |n|(/3o - ft+i + 1) ^" T'l-l'V 



i=l \ 1=1 



(2.11) 
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Similarly, we can write £p_2 ~ J2i<i<j<p ^p 'i where 



Using ()2.5p one can deduce as above that 

4:2^ (A^* ) = -2n,n,A^'- mod v1'^''^~^ 



which shows that 

'Cp.2(A") = -2 1 n,nj \ A" mod T'f'^V (2.12) 

\i<i<j<p / 

The proof follows easily from equations (|2.1ip and (|2.12p . □ 

3. MULTIVARIABLE RACAH POLYNOMIALS 

In this section we show that the operator Lp defined by ()2.6|) is self-adjoint 
with respect to the Racah inner product defined by Tratnik [27|- In order to 
unify the formulas in this section, besides the variables x\^xi, . . . ,Xp we also use 
xo = and Xp+i = iV. Thus, in view of equation (j2.ip . we can put Aq = and 

Ap+i = A^(A^ + /3p+i). 

3.1. The operator Lp in terms of the shift operators. The operator £p was 
defined in (|2.6p in terms of the forward and the backward difference operators 
{Aa;^ , Va;^ }j=iJ which was useful to prove the triangular structure of Lp. Below we 
give the explicit form of Lp in terms of the shift operators i?^, which allows us to 
show that Lp is self-adjoint with respect to the Racah inner product. 
For i e {0, 1, . . . ,p} and (j, k) G {0, 1}^ we define J5f as follows 

i..°^° = A. + A.,,+ (^^±lM±i-ll (3.1a) 
S"^^ = {xi+x + a:i + /3i+i)(xi+i -Xi+ Pi+i - f3i) (3.1b) 

■^i^'" = i^i+l - Xi){Xi+i +Xi+ Pi+i) (3.1c) 

bI'^ {xi+i +Xi+ Pi+i){xi+i +Xi+ Pi+i + 1). (3. Id) 
We extend the definition of Bj'''{z) for (j, k) G { — 1, 0, 1}^ by defining 

B,"i^'==/,(B^'^) forfc-0,1 (3.1e) 

Sf"' = /,+i(Bfi) for J = 0,1 (3.1f) 

Br'^-'^W,+x{Bl^')). (3.1g) 
Next, for i G {1, . . . ,p} we denote 

b1 = (2a;, + l)i2x, + - 1) ^ 4A, + pf ~ I (3.2a) 

51 = (2a;, + ft + l)(2a;, + ft) (3.2b) 

br' = h{b\). (3.2c) 
Finally, for G {-1, 0, 1}^ we put 

C. = 2^-l^-^l-k-| li/c=o^fe . (3.3) 
llfc=i Ofe 
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The main result in this subsection is the proposition below, which says that 
is the coefficient of i?^ for the operator Cp. 

Proposition 3.1. The operator Cp — Cp{x; /3; N) can be written as 

Cp{x;P;N)= Y: C.K-(\^i+ ^^°^'^f^'''^ ), (3-4) 
where Cu are given by (j3.1[) . (j3.2[) and (j3.3[) . 

Proof. From equations p.ip . (|3.2p and (|3.3p it follows that the operator £p{x; /3; A^) 
defined by (|3.4p is /-invariant. We prove the statement by induction on p. For p = 1 
formulas (|2.5p give 

. ^ (xi + /3i - /3o)(a:i + + ^2 + N){N - xi) 

(2a;i+/3i)(2xi+/?i + l) 
^ ^ xijxi + (3q){N - xi ~ I3i + f32){N + xi + (3i) 
(2a;i+/3i)(2xi+/?i-l) 
Using ([31]), ([33) and ((331) we get 

_ ^ (xi+/3i)(a;i+/3i-/3o)(jV-xi)(xi+/32+A^) ^ 

6} (2xi+/3i)(2:ri+/3i + l) 

_ . BgX'" _ + (/3o + - 1)/2)(A2 + Ai + + - l)/2) 
" " 4Ai+/3?-l ' 

where Ai — xi{xi + f3i) and A2 = N{N + (32). We need to check that 

which amounts to checking that 

This equality can be verified by a straightforward computation using the explicit 
formulas for A^-^, ^-ei and Cq. 

Let now p > 1 and assume that the statement is true for p — 1. Let us write Cp 
as follows 

r _r' {xp + Pp+i + N){N ~xp) 
" {2xp + Pp){2xp+(3p + l)^ 



-U ^ z. z.-!,^ A2 + (/?0 + l)(/32-l) 



(3 5) 

(A^ - a;p + /3p+i - /3p)(jV + + 

(2a;p + /3p)(2a;p + /Jp-l) ^ + ' 

where C , C" , C" are difference operators in the variables xi, a;2, . . . , Xp^\ with coef- 
ficients depending on xi, . . . , and the parameters /3o, . . . , /3p+i, N . Clearly, the 
operators CJ ^ C" , CJ" are uniquely determined from Cp. This implies that they are 
Ii, I2, ■ . ■ , Ip-i invariant and therefore they are characterized by the coefficients of 
= A'^\A'^l ■ ■ ■ AxlZ\, where x = {xi,X2, ■ . ■,Xp^i) and D ^ {i^i, 1^2, ■ ■ ■,i^p-i) £ 
{0, l}^'^^. The coefficient of A" in C is equal to —{xp + (3p — Po){xp + (ip) (it 
comes from the term A^^A^^ in Cp). Using equations (|2.5p - (|2.6p one can see that 
the coefficients of A^ for v G {0,1}^"^ \ {0} are the same as the coefficients of 
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the operator £p_i (xi , . . . , Xp-i; (3o, - ■ ■ , Pp-i, P'p\ N') where P'^ — 2{xp + (3p) + 1 and 
N' — —Xp — Pp. Thus we deduce that 

C =Cp-i{xi, . . . , Xp-.i;f3o, ■ ■ ■ ,/3p-i, 2(a;p + f3p) + 1; ~Xp ~ Pp) 
- (xp + Pp~ Po)ixp + Pp) 

iye{-i,o,i}p-^ 

where C'l are computed from ()3.H) - ()3.3I1 for the operator Cp-i with parameters 
given in ()3.6a|) . Notice that the last term on the right-hand side of p.4p for £p_i 
cancels with {xp + /3p — Po)ixp + Pp) in (I3.6ap . For C" we have 

For the last operator C" we obtain 

= £p_i(xi, . . .,Xp-i;Po, . . . ,Pp-i,Pp+i;N) 



-1 



where C"' are computed from p.ip - (j3.3p for the operator Cp^i with parameters 
given in (|3.6cp . 

Notice that the last term in (|3.6cp gives the last term in p.4p . Thus, it remains 
to show that we get the stated formulas for the coefficients C„ in p.4p . using the 
decomposition (|3.5p . Since the operator Cp is /-invariant, it is enough to prove the 
formulas for Ci, when v has nonnegative coordinates. We have two possible cases 
depending on whether Vp = 1 ov Vp = 0. 

3.1.1. Case 1: Vp ^ 1. Write v = {v' , 1) with v' e {0, 1}^"^ From (1^3]) it is clear 
that appears only in the first term on the right-hand side and we have 

_ {Xp + Pp+,+N){N-Xp) _ CI, , , ^ 

- {2xp + Pp){2xp + Pp + 1)- bl + + ^^^^ 

Notice that the factors (Sf in formula JSS]) for C'^, are the same as the 
factors in formula (|3.3p for for i = 0, 1, . . . ,p — 2. This combined with 

the last formula shows that in order to complete the proof we need to check that 

Again we have two possibilities: 

If Vp-i = then from the parameters in formula (I3.6ap we get 

(s^)' = Ap-i + N'{N' + p;) + (^-^ + m-^) 

^ Xp-i{Xp-l + Pp-l) - {Xp + Pp){Xp + Pp + l) + iPp-1 + l)(Xp -I- Pp) 

On the other hand from (|3.1bp we have 

Bp-l = {Xp + Xp^l + Pp){Xp - Xp-l + Pp- Pp-l). 

From the last two formulas it is easy to see that {B'^'\)' — —B'^'^. 
If Vp-i = 1 then 

(S^fi)' = {N'-Xp-i){N'+Xp-i+P'p) = -{xp+Xp-i+Pp)ixp+Pp+l+Xp-i) = -B^l 
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completing the proof in the case Vp = 1. 

3.1.2. Case 2: Vp = 0. Let us write again v = {v' ,Q) with v' e {0,1}^"^ Then 
from p.Sp we deduce 

^ {xp + Pp+, + N)iN - xp) 

" {2xp + f3p){2xp+(3p + l) 

^„ iN-Xp + Pp+i -l3p)iN + Xp + I3p) 

{2xp+(5p){2xp + f}p-l) ^^"'^ 

We need to check formula (|3.3p for C^. Again the factors B'^^''^^*^ for i = 0, 1, . . . ,p— 
2 and the denominator factors b^^ for i = l,2,...,p — 1 are common for C^, CI,, 
C", and C"l. Thus we need to verify that 

2i3;ri"°i?°'" _ , ^.^_„o,, {xp + Pp+i + N){N - xp) 



{2xp + f3p + l){2xp + pp - I) ' ' (2xp + /?p)(2a;p + /3p + l) 

{N -Xp + iip+i~(ip){N + Xp^(ip) 
^ ' {2xp + f3p){2xp + (3p-l) > ■ 

Using the explicit formulas (|3.ip and considering separately the two possible cases 
fp-i = and lyp-i = 1 one can check that the above equality holds, thus completing 
the proof. □ 



3.2. Racah inner product and polynomials. For a positive integer N we define 



an inner product on the space by 



if, 9)= E f{^)9{x)p{x), (3.7) 

where /, g € , Vm is the set 

Vat = {a; G Nq : = xo < a;i < a;2 < • • • < a;p < Xp+i ~ N} 
and the weight p{x) is given by 

Remark 3.2. If we put (3o = ai — r] — 1, fik = zZj=i fo'^ ^ ~ li 2, . . . ,p + 1, 
7 = — A^ — 1, then one can check that the weight defined in (|3.8p differs by a constant 
factor (independent of x) from the weight used by Tratnik, see [571 formula (2.3) 
on page 2338]. 

In [27j one finds an exphcit orthogonal basis in ■ If we denote by r„(a, /9, 7, 8] x) 
the one dimensional Racah polynomials 

rn{ct, f3, 7: S; x) ={a + !)„(/? + S+ 1)„ (7 + 1)„ 

(3.9) 



X 4F3 



— n, n + a + P + 1, -x, x + ^ + 6 + l_^ 
a + l,/3 + (5+ 1,7+ 1 
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then an orthogonal basis in Vj^ is given by the polynomials 

p 

Rp{n; x; /?; iV) = [] r„, {2N^-^ + /J, _ /?„ _ i, p^^^ - - 1, N^-^ - x^+i - 1, 
^1^"^ +I3k + xk+i; -N^-^ + xk). 

(3.10) 

Here /? = (/3o, /3i, • ■ • , /3p+i) are the parameters appearing in the weight (|3.8p . n = 
{ni,n2, ■ ■ ■ ,np) G Nq is such that |n| < N and iVj* = ni+n2 + - ■ ■ + nj with iV° = 0. 

Remark 3.3. Notice that the fc-th term in the product on the right-hand side in 
(IXTOll is 

p r-nfc,nfc + 2Ar^-i+/3fe+i-/3o-l,Aff~'-Xfe,iV^-i+/3fc+.Tfc " 

(3.11) 

From this formula it is easy to see that each r„^ is an /-invariant polynomial of x 
and therefore Rp{n\ x; /3; N) G V\. 

3.3. Self-adjointness of Cp. For a difference operator L = Si^es C^{x)E'^ (where 
5 is a finite subset of Z^') we can always assume that S is symmetric about the 
origin, by adding finitely many coefficients C^{x) identically equal to zero. In the 
next lemma we show that simple conditions relating Cv{x) and C-u{x) for every 
1/ ^ and certain boundary conditions allow us to check whether the operator L is 
self-adjoint with respect to the inner product (|3.7p . 

Lemma 3.4. // an operator L — X^z^es Cu{x)E'^ satisfies the following conditions 
for every ^ v E S 

(i) p{x)Cu{x) — p{x + v)C-y{x + v) when x,x + v ^ Vn ; 

(ii) Cv(x) = when x G Vn hut x + v ^ Vn 

then L is self-adjoint with respect to the inner product (13. 7p . 

Proof. For 7^ e 5 let Ly = Cy{x)E'' + C-y{x)E-''. Then 

{Luf,g)^ X! Cv^^)fi^ + ^)9{x)p{x)+ ^ C-^,{x)f{x ~ v)g{x)p{x) 

which becomes, using condition (ii), 

= X! C^{x)f{x + v)g{x)p{x) 

+ ^ C^^{x)f{x - v)g{x)p{x), 

replacing a; by a; — in the first sum and a; by a; + in the second 
= X! C^{x - v)f{x)g{x - v)p{x - v) 

x : x — u^ V/v , X G Vn 

+ ^ C_i.(x + iy)f{x)g{x + + i^), 

x+z^G V/v ,2;^ Viv 
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then using condition (i) 

+ X! C^{x)f{x)g{x + i^)p{x), 
again using condition (ii) 

xEVn x^Vn 

= {f,L,g). 

The proof follows immediately by writing L as a sum of L^^s. □ 

Lemma 3.5. For N N, the operator Cp{x] f3] N) is self-adjoint with respect to 
the inner product defined by p.7p - p.8p . 



Proof. Using Proposition 13.11 we show that the conditions in Lemma 13.41 are satis- 
fied. Let us write p{x) in ()3.8p as 

k=0 k=l 

where 

//N T{Pk+i- l3k + Xk+i-Xk)T{/3k+i+Xk+i+Xk) /o1o^ 

PkW = 7 ; ^ (3.12a) 

{xk+i-Xk)lT{Pk + l + Xk+i+Xk) 

pI{x) = + 2xk. (3.12b) 

Then condition (i) in Lemma 13.41 will follow if we show that for x,x -\- u G Vat we 
have 

p'ki-) ' B---^-^ix+.) ^'-''"^ 

and 

p'^x + iy) _ &^-'°(^ + ^) 

p'^i^) Kix) ■ 

Notice that equations (|3.13p for v and —v are essentially the same. Thus it is 
enough to check (jXTOal) for {vk,Vk+i) = (0,0), (1,0), (0,1), (1,1), (1,-1), while 
(|3.13b[) needs to be verified for j/fe = 0, 1. 

Let us start with (j3.13bp . li Vk = Q then clearly both sides of (j3.13bp are equal 
to 1. li Vk ~ 1 then from ()3.2|) we see that b\. — {2xk + /3k){2xk + Pk + 1) and 
b'l^^ = {2xk + (3k){2xk + — 1) which combined with (|3.12b|) shows that both sides 
of (|3.13bp are equal to {2xk + Pk + 2)/{2xk + Pk)- 

Next we verify equation ()3.13a|) . 
Case 1: Vk — Vk+i = 0. Clearly both sides of (|3.13ap are equal to 1. 
Case 2: Uk = 1, Vk+i = 0. From (13. 1|) we see that 



(3.13b) 



= {xk+i - Xk){xk+i +Xk+ Pk+i) (3.14a) 
Sfci^" = 4(B^'") = [xk+i + xk + Pk){xk+i ~xk + fik+i - Pk). (3.14b) 
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Using the well known property of the gamma function 

r(a;) = {x — n)„r(a; — n), (3.15) 

one can easily verify that the left-hand side of p.l3ap gives the same ratio as 
Bl^\x)/B-''\x + ,y). 

Case 3: = 0, i^fc+i = 1. In this case we have 

= [xk+i +Xk+ I3k+i){xk+i -Xk+ (ik+i - Pk) (3.16a) 

= Ik+i{B]:°) = [xk+i - Xk)[xk+i +Xk+ Pk), (3.16b) 

and applying p.lSp to the left-hand side of (|3.13ap one can check that it is equal 

to Bl-\x)/Bl -\x + I.). 

Case 4- Vk = Vk+i = 1- This time we need 

bI'^ = (xfe+i + Xk + /3fe+i)(xfe+i + Xk + (ik+i + 1) (3.17a) 

B^^^-^ = Ik+i{Ik{Bk')) = {xk+i + Xk + Pk)(xk+i + Xk + Pk - 1), (3.17b) 

and the verification of (|3.13a|) goes along the same lines. 
Case 5: Vk = 1, Vk+i = ^1- In this last case we use 

^k^^ = h+i{Bl'^) = {xk+i - Xk){xk+i -Xk-l) (3.18a) 

Sfe ^'^ = Ik{Bl'^) = (xfe+i -Xk+ Pk+i - Pk)ixk+i - Xk+ Pk+i - Pk + 1), 

(3.18b) 

and the verification of p.l3ap can be done as in the previous cases. 

Finally, we need to check condition (ii) in Lemma [3.41 Let ^ v E { — 1,0, 1}'' 
and let x G Vat such that x+i' ^ Vn- Since x e Vn it follows that Xi < Xi+i for every 
I G {0, 1, . . . ,p}. The condition x + i^ ^ Vn means that for some fc S {0, 1, ... ,p} we 
have Xk + Vk > Xk+i + Vk+i- Here again we use the convention that vq — Vp+i — 0. 
Since Xk < Xk+i the last inequality implies that Vk > Vk+i- Thus (vkji^k+i) must 
be one of the pairs: (0, —1), (1, 0), (1, —1). 

If Vk = 0, Vk+i = ^1 then we must have that Xk = Xk+i- Since C^{x) contains 
the factor which is when Xk ~ Xk+i (see (j3.16bl) '). we conclude that Cu{x) — 

0. 

If Vk — 1, Vk+i = then again we must have that Xk — Xk+i- This time Cu{x) 
contains the factor B]^'^ which is when Xk = Xk+i (see p.l4ap ') and therefore 
C^x) = 0. 

Finally, if Vk = 1, Vk+i = —1 then Xk — Xk+i or 1 + Xk — Xk+i- In both cases, 
B]:"^ = (see p.lSap ') and therefore C^{x) =0. □ 

3.4. Admissibility of Cp and the commutative algebra Ax- As an immediate 
corollary of Proposition 12.41 and Lemma [3.51 we obtain the following theorem. 

Theorem 3.6. Let iV e No and let {Q(n; A(x); /3; iV) : \n\ < N] he a family of 
polynomials that span such that Q{n; X{x); /3; N) is a polynomial of total degree 
\n\ which is orthogonal to polynomials of degree at most |n| — 1 with respect to the 
inner product p.7p - p.8p . Then Q{n] X{x); P; N) are eigenfunctions of the operator 
Cp{x\p\N) with eigenvalue — |n|(|n| — 1 -|- /3p+i — /3o). 

Remark 3.7. Following [3D] we can consider the ideal 

3(V/v) = {qeVx: q{x) = for every x G Vw}- 
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If we denote 

cr„(x) = {-Xi)m{-X2 + ni)„2 • • • {-Xp + ni + n2 -\ h np_i)„^, 

then it is easy to see that cr„(x) = when \n\ > + 1 and x G V^v- This shows 
that the /-invariant set 

&N = {(-a;i)ni(a;i + Pl)nA-X2 + ni)n^{x2 + ^2 + ^^l)n2 • • • 

X {-Xp + Til H h np^i)n^{xp + /3p + ni H h ?ip-i)np : |7i| = + 1} 

is contained in 3(V/v) and therefore for every q E Vx there exists q G such that 
q — q E 3(V/v). If n G Ng is such that |n| < N then the polynomial Rp{n; x; /3; N) 
defined by p.lOp has a positive norm which means that it cannot be identically 
equal to zero on Vn- Moreover, since {Rp{n;x; (3; N) : \n\ < N} are mutually 
orthogonal , we see that there is no polynomial (other then 0) of degree at most 
belonging to 3{Vn)- Thus the ideal 3{Vn) is generated by the set &n which allows 
us to identify the factor space V\/3{Vn) with ■ 

Theorem 13.61 shows that the operator Cp is admissible on Vn because for ev- 
ery fc g {0, 1, . . . , N} the equation CpU — ~k{k — 1 + Pp+i — Mu has C^^^Y^) ~ 
dim('P^'/'P^~^) linearly independent solutions in and it has no nontrivial solu- 
tions in the space V\~^ , see [16 , Definition 3.3]. 

Next we focus on the polynomials Rp{n; x; /3; N) defined by p.lOp for arbitrary 
(/?, N) E RP+'^, and we construct a commutative algebra generated by p difference 
operators which are diagonalized by Rp(n; x; (3; N). Before we state the main theo- 
rem, we formulate a technical lemma, where we prove the intuitively clear statement 
that a difference operator vanishing on V\ must be the zero operator. 

Lemma 3.8. If L = "^^^g Li'{x)E^ is a difference operator in W, with rational 
coefficients, such that L{q) = for every q E Vx, then L is the zero operator, i.e. 
Li,{x) = for every v E S. 

Proof. Without any restriction we can assume that 5 C Ng. Indeed if we fix ^ G Nq, 
then L is the zero operator if and only if i?^ • L is the zero operator. Thus, it is 
enough to prove the statement for the operator i?^ • L and choosing the coordinates 
of I large enough, this operator has only nonnegative powers of Er^. . 

Next, we can replace Ex^ by /\xj + 1, which shows that every operator L = 
X^i/es ^'^{^)-^x with 5* C Nq can be uniquely written as 



where M G N is large enough so that S' C 5a/ = {i^ G Ng : \v\ < M} and we have 
set L'^{x) = for J/ ^ S". 

Now we use the fact that L(A™) = for all m G Sm- We obtain a system of 
linear equations for Ll{x) with determinant 




\u\<M 



D = det 

u,m^SM 




We prove below that 



D = {2PXlX2 • ■ • 



M + p\ 

p+i ) + polynomial in x of total degree < p 
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which shows that D is a nonzero polynomial and therefore all coefScients Ly{x) 
must be equal to 0. 
Notice that 

^2m-v ^ polynomial in x of total degree < 2\m\ — \v\ 



I 



Moreover, 



where 











det 


= det 











M 



j=i 



M - j +p - 1\ (M + p 



p-l 



p+l 



Thus, to prove p.l9p it remains to show that 



det 



2m 

V 



(3.20) 



Let us denote C(M,p+ 1) = {a e Ng+^ : |a| = M}. Then by [F, Theorem 5, p. 357] 
for every z,l & we have 

(M + P\ 



det 

a,Q'eC(M,p+l) 



z + la 
a' 



n ..>o 

ci+e2 + ---+Cp+l<M 



where la denotes (Ziai, ^20^2, • ■ • , ^p+iQ^p+i)- If we apply this identity with a — 

(m, M — |m|), a' = {v, M — zi = Z2 = ■ ■ ■ = Zp = 0, li = I2 = ■ ■ ■ = Ip = 2, 
Ip+i — we obtain 



det 



2m\ / Zp+i 



n 



ei>0 ^ 



One can check now that 



det 



2to\ / Zp+i 
V ) \M ~ \v\ 



M 

n 



^p+i 



det 



2m 

V 



(3.21) 



(3.22) 



and 



n 

fc=i 



Zp+lV ^ _ ni=0 -J)^ " 



n e>>0 (2^p+l - Cp+l) 

£l+e2H hep+l<M 



(3.23) 



Equation p.20p follows immediately from p.2ip , p.22p and (|3.23p , thus completing 
the proof. □ 



From now on, we use the following convention: 
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• when ^ No we consider the polynomials Rp{n; x; /3; N) for all n e Nq; 

• when e No we consider the polynomials Rp{n;x; /3; N) for n G Nq such 
that \n\ < N. 

Theorem 3.9. For (/3, N) e Rp+^ and for j e {l,2,...,p} we define 

-^i = ■ ■ ■,Xj]l3o, . . . ,/3j+i;xj+i) (3.24a) 

Aij(n) = -(rii + n2 + ■ ■ ■ + nj){ni + n2 + • ■ • + — 1 + — (3o). (3.24b) 
Then 

Rp{n; x; P; N) = fij{n)Rp{n; x; (3; N) (3.25) 
for every j 6 {1,2, ... ,p} and the operators £J commute with each other, i.e. 
Ax — R-[£,f,£,2, ■ ■ ■ 1 ^p] is a commutative subalgebra ofT>x. 

Proof. If G No equation ()3.25|) follows from Theorem 13.61 Indeed, the product 
of the first j terms on the right-hand side of (|3.10p is precisely 

Rj{ni, . . .,nj;xi, . . .,Xj;po, . . . , Pj+i; Xj+i), 

while the remaining terms do not depend on the variables xi, . . . ,Xj. 

Notice next that Rp{n;x;(3;N) is a polynomial in N. Thus, if we fix n G Ng, 
then both sides of ()3.25p are polynomials in N . Since ()3.25p holds for every TV S Nq 
such that > \n\ we conclude that it must be true for every A^ G R. 

Finally, from p.25p it follows that for z, j e {1, 2, . . . ,p} 

£J]i?p(n; x; (3; N) = 0, for aU n e Ng. 

Thus, if A^ ^ No then [Z^ ,Z^]q = for every q ^ P\ which combined with Lemma 
13.81 shows that [£f,£|] = 0. Since the coefhcients of the operator [£f,£|] are 
polynomials of A^, it follows that if £f and £J commute for A^ ^ Nq they will also 
commute for A^ e Nq which completes the proof. □ 

Remark 3.10. From (|2.5p it follows that the operator Cp defined in (|2.6p is in- 
variant under the change A^ —N — (3p+i. This shows that all operators £J 
for j = 1,2, ... ,p, given in (|3.24ap are /-invariant, i.e. the whole algebra Ax is 
/-invariant. 

4. BISPECTRALITY 

4.1. Duality. Let x = {xi, X2, . . . , Xp) G W, n = {ni,n2, . . ■ ,np) G W, (3 = 
(/3o, I3i,..., (3p+i) G and let us define dual variables x G R^, n G and dual 
parameters (3 G R^"*"^ by 



X, - 7Vf+i-' + Pp+2-^ - /3o + iV - 1 for z = 1, . . . ,p, (4.1a) 

ni^ip + i3p + N, (4.1b) 

«j = Xp+i-i - Xp+2-t + i3p+i-i - Pp+2-i for i ^ 2, . . . ,p, (4.1c) 

Po = Po (4. Id) 

P^=Po- Pp+2-^ -2N + 1, for z = 1, . . . ,p + 1 (4.1e) 
with the usual convention that N( — hi + n2 + ■ ■ ■ + fij . 



Lemma 4.1. The mapping f : {x,h,P) —>■ {x,n,P) given by (|4.ip is a bijection 
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Proof. From equations (j4.1dj) - (|4.1e|) one can easily deduce that mapping 

{Po, Pi, ■ ■ ■ , Pp+i) — > (/3o, Pi, ■ ■ ■ , Pp+i) 

is one-to-one, and that the inverse is given by the same formulas. From (|4.1bp and 
(|4.1cp we find 

Ni^ip+i^i+Pp+i^i+N, ior i = l,2,...,p, (4.2) 

which shows that x can be written in terms of n and /3 by a formula analogous to 
(|4.1ap . From (|4.1ap with i ~ p wc obtain fii — Xp — N — P2 — Po + ^ = Xp + N + (3p. 
Finally, subtracting (14. lap for i and i+1 we get Xi — Xi^i — np-f-i^i+Pp+2-i — Pp+i-i 
which leads to a formula analogous to (|4.1cp for n^. □ 

Lemma 4.2. For k ^ l,p the polynomial r„^ in equation (j3.1ip can be rewritten 



r„, = (1 - iVf - Xk - /3fc)„fc (Pk+i - Pk)n, i-N^ - /3fc /3o + 1 - Xk)n^ x 

p \-nk, -Xk+i - Xk- Pk,Xk+i - Xk + Pk+i - /5fc, 1 - 2iVf - Uk ~ Pk + Po . 
^ U l~N^~Xk-Pk, Pk+i - Pk, -N^ -Pk+Po + l-Xk 

(4.3a) 

and for fc = 2, 3, . . . ,p — 1 we can write r„j. as 

Tuk = (1 - - Xk - Pk)nk [N'I^^ + Xk+1 + Pk+1 - Pa)nk (xk+l - Xk - Uk + l)nk 



-nk,Xk+i -Xk+ Pk+i - Pk, ^ " Po- Xk, 1 - TVf + Xk+i . ^ 
1 ~ - Xk ~ Pk,N^'^ + Xk+i + Pk+i ~ Po, Xk+i -Xk-Uk + l' 



(4.3b) 

Proof. First we iterate the identity of Whipple (see [21 p. 56]) connecting the ter- 
minating Saalschiitzian 4^F^ 



(M)n(w)„(w)„ 4F3 



-n,x,y,z ^ 



U, V, w 

= {1 - V - Z - n)n{l - W + Z - n)n{u)n 4-^3 



— n, u ~ x,u — y, z 
1 — v + z~n,l~w + z~n,u^ 



(4.4) 

Using this transformation again on the right-hand side of (|4.4p with z and u replaced 
by u — a; and I — v + z — n, respectively, yields 



{u)n{v)n{w)n 4-^3 



-n,x,y, z 



(4.5) 



X4F3 



1 



U, V, w 

(1 - a; - n)„(l - v + y - n)„(l - v + z - n),, 

-n, w — x,u — x,l — V — n 
I — x — n,l — v + y — n,l — v + z — n' 

Applying g31) with x ^ N^^'^ + Pk + Xk, y = Uk + 2N^-^ + Pk+i - Po - I, 
z = N^-^ -Xk,v = 2N^-'^ +f3k-po,u^ N^-^ + Pk+i + Xk+i, w = N^-^ - Xk+i 
we find for k — l,p equation (|4.3ap . 

When fc = 2, 3, . . . ,p - 1 we choose x = N^^'^ + Pk + Xk, y = Uk + 2N^^'^ + 
Pk+i - Po-1, z = N^-^ - Xk, V = N^-^ - xk+i, u - 2Nt^ + - Po, w = 
N^~^ + Pk+i + Xk+i and we obtain equation (|4.3bp . □ 

Using (|4.ip - (|4.2p one can deduce the following lemma. 
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Lemma 4.3. If {x,n,j3) and (x,n,(3) are related by (|4.ip and ifn^N^ then 
-ni, -X2 - xi - Pi,X2 - xi + /32 - 1 - ni - /3i + /Sq 



4F3 



4-F3 



1 - ni - xi - /3i , /32 - A , -ni - (3i + (Sq + I - xi 

'-N - - 1 - 27Vf -np-f3p+ /3o, -Up, N + (3p+i - Xp - Pp , ^ 
1 - iVf - £p - /3p+i - /3p, -TVf - ip - /3p + /3o + 1 



(4.6a) 



and for fc = 2, 3, . . . ,p — 1 we have 



4^3 



-nk,Xk+i - Xk + Pk+i - Pk.N'l ^ - /3o - Xfc, 1 - iVf + a;fe+i ^ 
1 - -Xk- Pk,N^~^ + Xk+i + Pk+i - Po, Xk+i -Xk-Uk + l' 



4F3 



Xp+2-k — Sp+l-fc + Pp+2-k — Pp+l-k, —flp+l-k, 



1-N: 



p+l-k 



ip+l-k — Pp+l-k,Nf + Xp+2_fc + (ip+2-k — Po, 



(4.6b) 



1 - ' + Xp+2-k,Nf ' - /3o - Sp+i-fc . 

Xp+2-k — ip+l-fc — ^p+l-fe + 1 

For n G Nq let us normalize Racah polynomials as follows 



Rp{n;x;/3;N) = 



Rp{n;x;/3;N) 



(4.7) 



(-7V)|„|(~iV - Po)\n\ nLi(A+i - Pk)r 
The main result in this subsection is the following theorem. 

Theorem 4.4. If {x,n,f3) and {x,h,/3) are related by (|4.ip and ifn,h £ Nq then 

Rp{n;x;P;N) ^ Rp{n;i;(3;N). (4.8) 

Proof. Let us write Rp(n]X](3]N) and Rp{h;x;P;N) using the representation of 
r„^, given in Lemma lL2l Notice that according to Lemma l473l the 4F3 factor in r^ 
is equal to the 4-^3 factor in r^^ (and therefore by Lemma l4.1l the 4_F'3 factor in r„p 
is equal to the 4F3 factor in r^j), and for = 2, 3, . . . ,p — 1 the 4F3 factor in r„j. 
becomes the 4^3 factor in r^^^j ,,. 

Thus if we compute the ratio Rp{n; x; f3; N)/ Rp{n;x; f3; N) all 4F3's will cancel 
and we get 

Rp{n;x;(3; N) ^ 
Rp{n;x;f3;N) 

TT (1 - -Xk- Pk)n, if3k+l - f3k)n, {-N^ - + /jp + 1 - Xk)n, 

kLi,p (1 - ^1^ -^k- h)n, 0k+l - Pk)n. {-Nl - + /3o + 1 - ik)n, 

p-1 



n 



(1 - jVf - Xk- /3fc)»,(^f + Xk+i + Pk+i - (3o)nAxk+i - Xk -nk + 1), 



k=2 (1 ~ ^1 -^k- Pk)hk{N^ + Xk+1 + Pk+i - Po)hk{xk+i -Xk-nk + l)ftj, 

(4.9) 
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Using (|4.ip - ()4.2p we can eliminate x and x on the right-hand side in (|4.9p and 
we can write it as the product T1T2T3, where 



^ IatP at »\ Vl 



Using 



for the numerator of Ti and an analogous formula for the denominator we see that 

^ _ {-N)ni ^ (~iV)|„| 
' (-iV)^P (-7V)|^I ■ 

For T2 we write 



p+l—k)nk 



Using that /3fe+i — Pk = Pp+2-k — /3p+i-fe and rearranging the terms one can check 
that 

^ (/3fe+l — /3fe)nfc 



Finally, a similar argument shows that 

{-N ~ /3o)jvP _ (~jV-/3o)|„| 



(-7V-/3o)^P (-iV-/3o)|. 



T3 = 



The proof now follows immediately from (14. 9p and the formulas for Ti, T2 and 
T3. □ 

4.2. The dual algebra An- We denote by T)^--^ the associative subalgebra of 
"Dx of difference operators with coefficients depending rationally on the parameters 
(/3, N). Clearly, the commutative algebra Ax defined in Theorem l3.9l is contained in 
T)^'^ . Similarly, we denote by T>!^'^ the associative algebra of difference operators 
in the variables n — (rii, ^2, • . • , with coefficients depending rationally on n and 
the parameters {P,N). 

Replacing i by p + 1 — fc in (|4.1cp we see that 

rip+i-k = Xk ~ Xk+i + (3k - Pk+i- 

From this equation it follows that a forward shift in the variable Xk will correspond 
to a forward shift in rip+i_fc and a backward shift in np+2-k for k — 2,3, ... ,p and 
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to a forward shift in rip when k = I. Thus, in view of the duahty estabhshed in 
Theorem 14.41 we define a function b as follows 

b{N) = N (4.10a) 

b(/3o) = l3o (4.10b) 

HPk) = /3o - /3p+2-fe - 2A^ + 1 for fc = 1, 2, . . . ,p + 1 (4.10c) 

b{xk) = ni + 712 + • • • + np+i^k + Pp+2-k -I3q + N -I (4.10d) 

b(i?,J = for fc = 1,2, . . . ,p (4.10e) 

with the convention that -Enp+i is the identity operator. 

Lemma 4.5. The mapping ()4.10|) extends to an isomorphism b : T)^'^ — T)^'^ . 
In particular, the operators £)! defined by 

£l = bi£t) for 1,2,..., p (4.11) 

commute with each other and therefore 

A„-b(A)-R[£?,£^,...,£;] (4.12) 

is a commutative subalgebra ofD^'^. 

Proof. Using (|4.10p one can check that 

b{E,,)-b{f{x)) = b{f{x + ek))b{E,J, 

holds for every fc = 1, 2, . . . ,p, which shows that b is a well defined homomorphism 
from T)^'^ to T>!^'^ . The fact that b is one-to-one and onto follows easily. □ 

For j — 1,2, ... ,p we denote 

Kj{x) = b-^i^i-jin)) = -{xp+i-j - N){xp+i-j + Pp+i-j + N) 
= -K+i-j{xp+i-j) + N{N + (3p+i^j). 
The main result in the paper is the following theorem. 

Theorem 4.6. The polynomials Rp(n;x;(3;N) defined by equations p.lO|) and 
()4.7p diagonalize the algebras Ax and An- More precisely, the following spectral 
equations hold 

£,^j Rp{n; x; P; N) ^ fij{n)Rp{n; x; l3; N) (4.14a) 
£]Rp{n;x;P;N) = Kj{x)Rp{n;x; (3; N), (4.14b) 
for j — 1,2, ... ,p where £|, /ij are given by p.24p and Hj are given in (|4.1ip . 

Remark 4.7 (Boundary conditions). Since £" contains backward shift operators, 
and since Rp{n; x; (3; N) is defined only for n G Nq it is natural to ask what happens 
when £" produces a term with a negative Ui for some i. From (|4.10ep it follows 
that b{E'^) with v G { — 1,0,1 will contain a negative power of En o^^c of 

the following two cases: 

Case 1: Vk = 1, Vk-i = 0. In this case, b(£'^) will contain E^^^^_^. Notice that 
the coefficient of E!^ is C^, which has the factor {xk — Xk~i + Pk — Pk-i) (in B'^'^^ - 
see formula p.l6ap ). Since b{xk — Xk-i + /3fc — Pk-i) — —"■p+2-fc we see that this 
term is when — 0. 

Case 2: Uk = 1, ^k-i = —1- This time biE'i) contains „ , . The coefficient 

\ jj / I tp_j_2 — k 
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C,y has the factor B^^'^ = {xk - Xk-i + Pk ~ I3k~i){xk - Xk-i + Pk - Pk-i + 1) 
(see (j3.18bp ). Since b^Bj^}'^) — npj^2-k{np+2-k ~ 1) we see that this coefficient is 
when = or 1. 

Proof of Theorem \4-6\ Equation (|4.14ap follows immediately from Theorem l3 . 9l and 
the fact that Rp(n; x; /3; N) and Rp(n; x; /3; N) differ by a factor independent of x. 
It remains to prove (j4.14bp . 

Consider first the case when N ^ Nq. Then Rp{n-, x; j3; N) are defined for all 
n S Ng. Let us fix n,/3,N. By Lemma [4.11 for every a; G there exist unique 
{x,n,j3) such that f(i,n, /3) = {x,n,P). If h £ Nq then we can use Theorem 14.41 
and replace Rp{n; x; (3; N) by Rp{n; x; /3; N). Equation (j4.14bp follows from the fact 
that the operator £" in the variables n with parameters (3, N coincides with the 
operator £J in the variables x with parameters /?, N. Now we can think as follows: 
fix n, p, N and write x in terms of the dual variables h, i.e. we put 

Xk = fii + n2 -\ h np+i-fc - /3fe - iV for fc = 1, 2, . . . ,p. 

Both sides of equation (|4.14b[) are polynomials in n of total degree at most |7i| + 1. 
Since (j4.14bp is true for every n S Nq, we conclude that it must be true for arbitrary 
ri G RP, or equivalently, for arbitrary x € R^. 

If G No we can obtain equation (|4.14bp by a limiting procedure. In this 
case, we consider the polynomials Rp{n]x-P;N) for n e Ng such that |n| < N 
and equation (j4.14b|) holds when \n\ < N. When \n\ = N all terms in (j4.14bj) 
have well defined limits, but the right-hand side will contain also polynomials of 
degree iV + 1. Indeed, from (|4.10ep it follows that b{E^) = E^^ where = or 
= ±1, and \v'\ = 1 if and only if vi = 1. In this case has xi + /3i as a factor 
(coming from B^ , see (|3.1b|) ). Since b{xi + = |n| — iV this term is equal to 
when \n\ — N. Notice that the polynomials Rp(n;x;P;N) are well defined even 
when |n| > iV-t- 1 (since all denominators in p.9p coming from the hypergeometric 
series cancel) . From (14.71) it follows that if G Nq and if n. A; G Nq are such 
that \n\ ^ N, |fc| = + 1, then limAr,^jv(|"-| - N')Rp{k; x; /?; N') is a well defined 
polynomial. Thus, each term on the left-hand side of (I4.14bp has a natural limit 
when TV G No. □ 

Remark 4.8. Using the isomorphism b in Lemma [4.5l we can also define involutions 
/" on the dual algebra V!^'^ by 

I^'^boIjob-\ for j l,2,...,p. 

Then the operators i2", i — 1,2, ... ,p will be invariant under the action of the 
involutions j = 1,2, ... ,p. 

In the one dimensional case, this involution played a crucial role in the construc- 
tion of orthogonal polynomials satisfying higher-order differential or g-difference 
equations, by applying the Darboux transformation to the second order difference 
operator £" corresponding to the Jacobi polynomials llj and the Askey- Wilson 
polynomials [T^. 

5. Examples of other bispectral families of orthogonal polynomials 

It is well known that the Racah and the Wilson polynomials are at the top of 
the Askey-scheme [19]. In this section we show that a change of the variables and 
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the parameters leads to bispectral algebras of difFerence operators for multivari- 
able Wilson polynomials. We consider also in detail several limiting cases related 
to multivariable Halm, Jacobi, Krawtchouk and Meixner polynomials, illustrating 
their bispectrality. 

5.1. Wilson polynomials. If we set a = a + b— 1, f3 — c + d— 1, 'j = a + d— 1, 
5 — a — d and x — —a — iy in (|3.9p we obtain the Wilson polynomials 



Wn{y; a, b, c, d) = {a + 6)„(a + c)„(a + d)„ 
-n, n + a + b + c + d 



X4-F3 



l,a + iy,a 
b, a + c, a + d 



ly _ 



If we make the change of variables 

= a — b 

Pk = 28^ + 2a 

^2 +2a + c + d 



Xk = -£2 - a-iyk 
N = -El -a-d, 



(5.1) 



(5.2a) 
(5.2b) 
(5.2c) 
(5.2d) 
(5.2e) 



where k = 1,2, . . . ,p, £2 — ^2 + ^3 + ■ ■ ■ + ^k, ^2 — 0, we see that the Racah poly- 
nomials Rp{n; x; (3; N) defined by (|3.10p transform into the multivariable Wilson 
polynomials Wp{n; y; a, b, c, d, e) defined by Tratnik , see [26l formula (2.1) on page 
2066]. When 

Re(a, b, c, d, €2, £3, . . . , ep) > 
these polynomials are orthogonal in MP with respect to the measure 

dp'aAc,dJy) = r(a + iyi)r{a - iyi)r{b + iyi)r{b - iyi)x 
p-i 



n 

fe=i 



r(efc +i + iyk+i + iyk)^{ek+i + iyk+i - iyk)^{ek+i - iyk+i + zyfc)r(efc+i - iyk+i - iyk) 

T{2tyk)T{-2iyk) 



r(c + iyp)T{c - iyp)T{d + iyp)T{d - iyp) 
T{2iyp)T{~2iyp) 

From (|5.2p it is clear that the shift E^f. corresponds to the shift operator Ey^ in 
the variables y defined by Ey^f{y) = fiy + isk)- Thus the bispectral properties 
established in Theorem l4.6l lead to bispectral commutative algebras Ay and An for 
Wp{n; y; a, b, c, d, e) after the change of variables (|5.2p . 

5.2. Hahn polynomials. Following Tratnik [57] we can obtain the multivariable 
Hahn polynomials of Karlin and McGregor [18] as follows. We introduce new 
parameters {71, 72, ... , 7p+i} by 

Pk=Po+li+k, for A: = 1,2, . . . ,p+ 1, (5.3a) 

where 7^ = 71 + 72 + ■ • • + 7*:, 7i = and new variables yi,y2, ■ ■ ■ ,yp hy 

Vk^Xk-Xk-i for fc= l,2,...,p, (5.3b) 

where xq — Q. Then 

- = 7fe+i + 1 for fc 1,2, . . . ,p. 
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Up to a factor independent of x, the weight p{x) in (13. 8p can be written as 

// N TT (/3fc+l - /?fc)a:fc + i-Xfc(/3fc+l)xfc+i+Xfc -pr /3fc + 2a;fc 

Thus, if we let /3o ^ oo we see that p'{x) approaches 

P. Ay) = n ^'7;"'?m'" = n (5-4) 

fc^i 2/fc! (A^ - \v\y- Vkl 

with the convention that Up+i = N — \y\. For a, 6, iV G K. and n G No we denote by 
hn{x; a; b; N) the one dimensional Hahn polynomials 



hn{x; a; b; N) ^ {a + l)„(-iV)„ 3F2 



—n,n + a + b+l,—x 
a + l,-N ' 



Then when /3o ^ 00 the polynomials Rp{n; x; f3; N) defined by p.lOp and (|4.7p 
become the multivariable Hahn polynomials 

(-1)1"! 



(-^)|n|nLi(7fe+i + l)n. 

(5.5) 

X n - Nfi^Nt' + 7i' + fc - 1; 7fc+i; - ^1'"'), 

fe=i 

orthogonal for iV G Nq with respect to the weight p^^N{y) given in (|5.4[) on the set 
y = {?/ G Nq : |j/| < A^}. Here we use again the convention = yi + y2 + ■ ■ ■ + Vk- 
Formula (|5.5|) above coincides (up to a normalization factor) with the formulas given 
by Karlin and McGregor, see [HI p. 277]. Since the weight is invariant under an 
arbitrary permutation of the labels (l,2,...,p+l), one can easily generate different 
orthogonal families, corresponding to different ways of applying the Gram-Schmidt 
process in the vector space of polynomials of total degree fc, for fc = 0, 1, 2, . . . , iV. 

Let us now focus on the difference operators {'CJ,'C"}j^i and the eigenvalues 
{/ij(ri), '«j(a;)}^^j^ in Theorem 14 . 61 when we change the parameters and the variables 
according to (|5.3p and then let (3q — > 00. First we change the parameters using 
(I5.3ap . We can write 

^= E ci^\P^,i,N){E:-i), 

where ci''\/3o,'y, N), B^'""''^^ {Po,j, N), etc. denote the coefficients and their com- 
ponents in the new parameters f3o,ji, ■ ■ ■ ,"fp+i, N . From the explicit formulas 
(I3ll])-(l3l2]) (see also ([XT4| - ([3?T8l) ) it is easy to see that 

^.^ B'^"'''+\f3o,-f,N) _ if = = 

Po^oc /3p 1 1 if z-'i = z^i+i = ±1 

hm IS unite when i^i ^ Vi+i, 

lim . 1. 



and 
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This combined with (13. 3|) shows that if for i = 0,1, ... ,j we have at least three 
pairs {vi,v.ij^i) such that ^ //j+i then \mifjg_,ooCy\l3Q,"f,N) — 0. Thus we 
need to consider only v ^ Q oi the form: vi = vi+i = ■ ■ ■ ^ Vk = il for some 
1 < / < A: < j and all other coordinates of v are equal to 0. Notice also that the 
operator ExiE^^^-^ ■ ■ ■ E^^ corresponds to the operator Ey^E~^^_^ in the variables 
yi,y2,. .. ,yp. Here we adopt the convention that Ey^^-^ is the identity operator. 
If we put m = k + 1 then after the change of variables and parameters (|5.3[) the 
operator £J reduces to the operator 

E iyi + ii + '^)ymiEy,E;^-i). (5.6) 

l<l=im<j + l 

This operator plays a crucial role in the work T8J. Another useful representation 
of £j is given below 

i+i 

= - E +^' + l)2/mA,, Vj„„ + J2 ym{Yl+^ - y^)Ay,^^y^ 

l<l^m<j+l m=l , , 

+ E + - ym)^y^ - E - lm+j)ymyy,„, 

with the convention that Aj^^^^ and Vy^^^ are zero operators, or equivalently the 
sums are up to p when j = p. 

Next we consider the difference operators in n. Let us denote by Li>'''' the coeffi- 
cients of the operator £" — £"(/3o, 7, N) defined by (|4.1ip in the variables {(3o,1, N), 
i.e. we can write 

£]{f3o,j,N)^ ^ L«(b(i?,^)-1), 
07^i/e{o,±i}3 

where 

4^') = b(C,.) = 2J-I''^I-I''"I . (5.8) 



On the right-hand side of the last formula, we have Bj, ' and b/!' which denote 
j^^^i^fc,'^fc+i-j i'ib^.''), respectively in the parameters /3o,7, A^. 

Lemma 5.1. The product 

3 ^Vk,Vk+i 



n 



k=l "k 

is independent of Pq and N, i.e. all terms in (jS.Sp except B^''^^ will not be affected 
by the limit /3o — > oo. Moreover, 

lim — ^ — = N + terms involving n and 7. (5.9) 

Proof. Notice first that fa(/3fc) and b{xk) are independent of Pq after the change of 
variables (|5.3p . Using also the fact that b(a;i+i —Xi), b(/3i+i — /3i), b{xi+i+Xi+f3i+i), 
b(a;i+i + Xi + Pi) are independent of A'^ for i > 1, we obtain the first assertion of 
the Lemma. Equation (|5.9p can be easily checked by considering the three possible 
cases vi = 0, 1, —1. □ 
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Based on the above Lemma we define 

1 



lim _£"(/3o,7,7V), forj = l,2,...,p, 



(5.10) 



where in the right-hand side the operator is written in terms of /3o and 7, using (j5.3p . 
Then as a corollary of Theorem 14.61 we obtain the following bispectral property of 
the multivariable Hahn polynomials. 

Theorem 5.2. The polynomials Hp(n; y; 7; TV) defined by (|5.5p satisfy the following 
spectral equations 

~S,)Hp{n- y- 7; N) = (n)i/p(n; y; 7; TV) (5.11a) 
£Ji?p(n;y;7;TV) = k,{y)Hp{n-y;r,N), (5.11b) 
/or J = 1, 2, . . . ,p where £| and £" are given by ()5.6p anrf ()5.10p . respectively, and 
jljin) = -(ni + ^2 H h nj)(ni + 7^2 H 1- + J + 7i H + 7i+i) 



= ^ - yi - 2/2 



Remark 5.3. It is interesting to notice that while the operators £^ are signifi- 
cantly simpler than £| (for instance, £| has p{p+ 1) terms with coefficients which 
are quadratic polynomials of y, while £p has 3^ terms with coefficients which are 
rational functions of x) the operators £" and £" have exactly the same number 
of components, and they both have rational coefficients of n. This phenomenon 
appears even when we take one more limit and we consider the Jacobi polynomials. 
This is the content of the next subsection. 

5.3. Jacobi polynomials. From the multivariable Hahn polynomials we can ob- 
tain multivariable Jacobi polynomials as follows. We introduce new variables 
zi,Z2,...,Zp by 

yk^Nzk, forfc = l,...,p. (5.12) 
Using the conventions in the previous subsection, we define Zp+i = 1 — |z| and 



Zi +Z2 H 1- 



huk in (|5.5p satisfies 
fim 

N~.oo TV"fc 



Zfe. Applying the substitution (|5.12p . we see that the fc-th term 



= (-Zf+i)"'=(2TV* 



X 2^1 



-rik-nk 



-r7iN 

rfc-l 



-7i 



k+i 



2iVf-^+7f + A: 'Zi^+i 
Let us denote by P„(x;a,/3) the usual Jacobi polynomials 



(5.13) 



P„(x;a,/3) 



(a -f 1)„ 



'Fx 



-n, n + a + (3 + 1 
a+1 



Using (|5.13p we see that if we make the substitution (|5.12p in (15. 5p and let — > cxo 
the multivariable Hahn polynomials become the multivariable Jacobi polynomials 



Jp(n;z;7) = (-1)1-1 n 



ukKz^ 



k=l 



ilk 



+ 1 



P 



1 



zi 



fc+ 



These polynomials are orthogonal on the simplex 
TP {z G : > for i = 1, 2, . . . ,p and \z\ 



(5.14) 



zi + ^2 + ■ ■ • + < 1} 
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with respect to the weight 

P-riz)^f[z^,'il-\z\r'^^K (5.15) 

k=l 

Notice that 

hm NAy, = hm NVy, = A. 

Thus, using (|5.7[) we see that if we make the change of variables (|5.12p and let 
N ^ oo, the operators £| become the differential operators 

q2 minO"+l,p) 

l</<m<minO + l,p) ^ 



min(^-|-l,p) 



(5.16) 



^ ((7,„+i)^r'-(7r'+j+i)z™) 

m—1 

When j — p we obtain the operator 



dzr, 



52 P ^2 

l<l<m<p m=l 



E ((7m + l)-(7r'+p+l)z™) 



m—1 



dZrr. 



This operator is independent of the way we apply the Gram-Schmidt process in 
each vector space of polynomials of total degree k, for k = 0, 1, 2, ... , i.e. if we 
construct orthogonal polynomials on with respect to weig ht (f5?T5l) they will 
be eigenfunctions of £p. The operator £p was derived at the beginning of the 
last century in the monograph [l] in the case p = 2, from the differential equations 
satisfied by the Lauricella functions. The operator for general p can also be obtained 
by similar techniques, see |17j . Other proofs use Dunkl's differential-difference 
operators, see 9, for details. 

We define difference operators in n by 

lim (5.17) 

By Lemma [5.11 these operators will have the same supports as the operators 
More precisely, using (|5.8p we can write 

05^i/6{0,±l}J 

where 



fc=l °fe 

Combining the above formulas with Theorem 15.21 we obtain the bispectrality of 
the multivariable Jacobi polynomials on the simplex with weight function given 
by (I51B . 
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Theorem 5.4. The polynomials Jp{n]Z]j) defined by (|5.14|) satisfy the following 
spectral equations 

£,jjp{n;z;^) = fij{n)Jp{n; z;j) (5.18a) 

£]Jp{n;z;j) = kj{z)Jp{n;y;j), (5.18b) 

for j = l,2,...,j3 where £| and £" are given by (|5.16p anc? (|5.17p . respectively, 
and 

AiC*^) = -(^1+^2 + h nj){ni + 7i2 H + j + 7i H + 7j+i) 

Kj(z) = 1 - Zl - Z2 Zp+l_j. 

5.4. Krawtchouk and Meixner polynomials. Multivariable Krawtchouk and 
Meixner polynomials can be obtained from the multivariable Hahn polynomials 
which is one way to construct the bispectral operators. On the other hand, 
they possess interesting duality properties, which combined with the explicit form 
of the admissible difference operators [16] , allows us to write simple formulas for the 
commutative algebras on both sides. Below we sketch the main ingredients of this 
construction, which parallels the theory developed for the Racah polynomials in the 
first sections. For interesting applications of Krawtchouk and Meixner polynomials 
to multi-dimensional linear growth birth and death processes see |22| . 

For n € Nq and parameters p = (pi,p2, • ■ ■ ,Pp) and N, we define multivariable 
Krawtchouk polynomials by 

Kp{n;x;p;N) = ^^—f[l^^ Ix,, '^^ ; N - \n\ + Nf - Xr'] , (5.19) 

where as before, for y Cz we set — yi + y2 + ■ ■ ■ + yj, with the convention 
Yi = 0, and kn are the one dimensional Krawtchouk polynomials 



knix;p;N) = i-N)n 2F1 



for n e No and x,p,N £ 



-n, —X 1 
-N ' p 

The multivariable polynomials defined by ()5.19p are eigenfunctions of the operator 

p 

Cp{x]P]N)^ p.XjA^^V^^. +^[p,(a;, -iV)A^^ + (l-p,;)x,V:,J (5.20) 

l<j#i<p i=l 

with eigenvalue fJ,p{n) ~ \n\. When N £ Nq and |p| < 1 these polynomials are 
orthogonal on {x e Nq : < N} with respect to the weight 

1 ^ 1 / pA: ^ 

^^""^^ iN-\x\y.Jl^.[T^\ 

Remark 5.5. If we put 

Pk = and N - -.s (5.21) 

\c\ - 1 

in (|5.19p we obtain the multivariable Meixner polynomials A4 (n; x; c; s) . When 
and s are positive with |c| < 1, these polynomials are orthogonal on Ng with respect 
to the weight 



PcA^) = (s)|.| [[ —f-. 



Xl.\ 

fc=l " 
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^ vp+i^/c.^ in; (5 22c) 



Thus all difference equations derived for the Krawtchouk polynomials will be true 
for the Meixner polynomials after changing the parameters according to (I5.21|) . 

Let us now introduce dual variables and parameters for the Krawtchouk poly- 
nomials by 

= fip+i-k (5.22a) 
rifc = Xp+i_fc (5.22b) 

Pp+i~fc(i - IpI) 

(l_^P+i-fe)(l_(pP-«) 

for k — 1, 2, . . . ,p. Then one can deduce the following analog of Lemma [4.11 and 
Theorem 14.41 In this case the duality follows immediately and we do not need any 
transformations for the 

Theorem 5.6. The mapping f : (i,n, p) — > (a:,n,p) given by (j5.22p is a bijection 
R^'f such that = f. Moreover, ifn,n € Nq then the polynomials defined 
by (j5.19p satisfy the following duality relation 

Kp{n;x;p;N)^Kp{n;i;p;N). (5.23) 

Notice next that if we consider in (|5.19p the product of the terms for k ~ j 
to p, then we obtain (up to an unessential factor) the multivariable Krawtchouk 
polynomial 

pj pj+1 pp . ^ x^^A 

Since the product of the terms for fc = 1, 2, . . . , j — 1 is independent of the variables 
Xj, Xj+i, . . . ,Xp we define for j = 1,2, . . . ,p the operators 

£^-r^, (x-x^^ X- Pj+i h -N-X^^A 

i.^ - L,p+i-j \Xj,Xj+i...,Xp, I. 

(5.24) 

Thus we obtain the commutative algebra 

Ax = M[£i , £2 J • ■ • ; -^p]- 

If we denote by V^'^ and V^'^ the associative algebras of operators in the variables 
X and n, respectively, with coefficients depending rationally on the parameters p, N, 
we define an isomorphism b : V^'^ T^n^ by 

b{xk) ^ Up+i^k, b{Ex^) = En^^^_^ k = \,2,...,p (5.25a) 

Pp+i-fc(i - IpI) 

(l_<pP+i-fc)(l_«pP-'=) 
b{N) = N. (5.25c) 
This leads to the dual commutative algebra 

An = ^2, ■■■,^p], where = b{£^). 



b{Pk) = fp+i-A,y iFM fc = l,2,...,p (5.25b) 
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Theorem 5.7. The polynomials Kp{n; x; p; N) defined by equations ()5.19p diago- 
nalize the algebras Ax and An ■ More precisely, the following spectral equations hold 

£,'jKp{n; x; p; N) = {uj + nj+i H h np)Kp{n; x; p; N) (5.26a) 

£,'JKp{n;x;p;N) = {xi + X2 H ^ Xp+i^j)Kp{n;x;p;N), (5.26b) 

for j = 1,2,..., p. 

Appendix A. Explicit formulas in dimension two 
A.l. Racah polynomials. From (|3.10p and (|4.7[) we have 
R2{n-x- (3-N) = '^"i ~ Po-l,P2- Pi~ 1, ~X2 - 1, /3i + X2; xi) 

( — A^)ni+n2 ~ Po)ni+n2{P2 — Pl)ni{Pz — ^2)712 

X r„2(2ni + 132 - Pq - 1, 133 - 132 - l,ni - N - l,ni + /?2 + N; -rii + X2). 
The operator £| can be written as follows 

£2= E c(2)(i?^-iij^-^-i) 

05^iye{-l,0,l}^ 



where the coefficients are given by 
C, 



(2) _ (Xi + Pi){xi l3o){x2 +Xi+ I32){X2 +Xi+(32 + 1){N ^ X2){N + X2+ Ps) 



(^■1) (2x1 + Pi){2xi +(3i + l)(2a;2 + /32)(2a;2 + /32 + 1) 

^(2) _ (xi + /3i)(xi +I3i- Po){x2 - xi){x2 +xi+ /32)(2A2 + 2A3 + {132 + l)(/33 - 1)) 
{2x1 + (3i){2xi + 1){2X2 +p2 + l){2x2 + f32 ~ I) 

M2) _ (2Ai + (/3o + l)(/3i - 1))(X2 +XI+ P2){X2 - .xi + /32 - /3i)(A ~ a;2)(Af + X2 + fe) 
(2x1 + /3i + l)(2a;i + /3i - l)(2x2 + /32)(2a;2 + 132 + 1) 

^(-1,1) - ^H<-^(l,l)i' "-^(l.-l) - -'2l'-(l,l)j' "-(-l.-l) - -'U-'2(<-(l,l)jj 
ri(2) _ T \ ^(2) _ J fM2) \ 

where Ai = xi(a;i +^1), A2 = X2{x2+l32), A3 = N{N + P^). 
Similarly we have 

&l^c['\Ex,-l) + d'l{E-'-l), 

where 

^(1) ^ (xi + /3i - /3o)(a:;i + f3i){x2 + xi + (32){x2 - Xi) ^(1) ^ ^ .^(1). 
1 (2x1 +/3i)(2xi +/3i + 1) ""'^ n 1 j- 

The corresponding eigenvalues are 

H2{n) = -(ni + n2)(ni + n2 - 1 + /53 - /^o) 
/ii(n) = -ni(ni - 1 + /32 - Pa)- 

Next we write the explicit formulas for the dual operators £2 ^^'^ ^i- First we 
write 



BISPECTRALITY OF MULTIVARIABLE RACAH-WILSON POLYNOMIALS 31 

where S = {(1, 0), (0, 1), (-1, 2), (1, -1), (-1, 1), (0, -1), (1, -2), (-1, 0)} and the 
coefficients are given by 

r,(2) ^ (m + na - N){ni + n2 - N - /3o)(2n.i + na + A3 - A) - l)(2ni + na + /?3 - /?o) 
(2ni + 2n2 + - Po - l)(2ni + 2n2 + /Js - Po) 

^ (2ni + /32 - /?o - l)(2ni + /32 - /3o) 



(2) _ (ni + 712 - jV)(ni + 712 - iV - /3o)(-n2 + P2 - PsK^rn + n2 + P3 - Po - 1)^2(71) 
(2ni + 2712 +P3-P0- l)(27ii + 2n2 + Ps - Po){2ni + P2 - Po){2ni +p2-po-2) 



' — 



^(2) ^ (7li + 7l2 - jV)(7ti + 7l2 - TV - /3o)(-7Z2 + P2 - Ps) 
(-1>2) (2711 + 2712 +P3-P0- l)(27ll + 2712 + Ps - Po) 

^ i~n2 +P2-P3- l)ni(-7ii +P0-P1 + I) 
{2ni + p2 - Po - l){2ni + p2 - Po - 2) 

^{2) ^ To(7i)(27ii + 772 + /33 - /3o - 1)^2(711 + P2 - Po - l)(ni +P2- Pi) 

(I'-i) (2711 + 2712 + Ps- Po){2ni + 2712 + P3-P0- 2)(2ni + P2 - Po - l)(2ni + /32 - /?o) 

^(2) ^ ro(7l)(-7t2 + /?2 - /?3)(277l + 7X2 + /32 - /3o " l)ni(-7ti + /gp - /3i + 1) 

(-1-1) (2ni + 2712 + /33 - /3o)(27ii + 2712 + /?3 - /3o - 2)(2ni + P2 - Po - l)(2ni + /32 - /3o - 2) 

(2) _ {ni + 772 + A3 - /3o + iV - l)(77i +n2+P3 + N~ l){2ni + n2 + P2 - Po - l)n2T2{n) 



D 



(O'-i) (2771 + 2772 + P3-P0- l)(277i + 2772 + P3 - Po - 2)(277i + P2 - Po){2ni +P2-P0-2) 
^(2) ^ (771 +n2 + p3-Po + N- l)(77i + 772 + /^s + " 1)^2(7^2 - l)(-ni - /?2 + /^o + l)(ni + /32 - Pi) 

(2771 + 2772 +P3-P0- l)(277i + 2772 + /33 - /3o - 2)(277i + P2 - Po - l)(277i + /32 - /3o) 

^(2) ^ (m + 772 + As - /3o + - 1)(tX1 +n2+P3+N- l)(277i + 772 + /^a - /3o - 1) 
(2771 + 2772 + /33 - /30 - l)(277l + 2772 + p3 ' Po ' 2) 
(2771 + 772 + /?2 - /3o - 2)771 (-771 + /3o - /3l + 1) 

(2ni + /32 - /3o - l)(2ni + /32 - /?o - 2) 
where we have set 

To(77) = b{2B^o'°) = -2M2(n) + {Po + l){Po - Ps) - 2N{N + ps) 
T2{n) = b(2B°'°) = -2/7i(n) + {po - pi){Po -p2 + 2). 
For £,1 we have 

= i)g;i)(i;„, - 1) + i)!o!_i)(i5-' - 1), 

where 

^(1) ^ (m + 772 - jV - /?o)(ni + 772 - A^)(2ni + 772 + /?3 - /^o - l)(-n2 + /32 - Ps) 
(2ni + 2772 + P3-P0- l)(2ni + 2n2 + /Ss - Po) 

and 

„(1) ^ (771 + 772 + /33 + ^ - l)(rtl +n2+P3-Po+N - 1)772(2771 + 772 + /32 - /3o - 1) 

(277i+2n2 + A3-/9o-l)(2ni + 27i2 + ^3-/3o-2) 
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Finally the eigenvalues are 

H2{x) = -\i{xi) + N{N + fh) 

= -\2{X2) + N{N + (52). 

A. 2. Jacobi polynomials. When p = 2 formula (|5.14p gives 

, , (-l)"i+"^ni!n2!(zi + Z2)"^ 
J2(n;z;7) = 



X P„, ( 1 71,72 ) Pn^Cl - 2(zi + Z2); 2rii +71 +72 + 1,73)- 



(72 + l)ni(73 + l) 

2zi 
2:1 + ^2' 
On the z side we have 

92 g2 ^2 

a d 

+ ((71 + 1) - (71 + 72 + 73 + 3)zi) — + ((72 + 1) - (71 + 72 + 73 + 3)2:2) 

^r.z n 5' 9' ^2 

£1 = -2ziZ2t; — ^1 h ziZ2Tr^ + 2122^7^ 

aziaz2 ozf az^ 

d d 

+ ((71 + l)-22 - (72 + 1)21) ^ + ((72 + 1)21 - (71 + 1)22) 

Next we write the explicit formulas for the dual operators £3 and We can 
write 

where S = {(1, 0), (0, 1), (-1, 2), (1, -1), (-1, 1), (0, -1), (1, -2), (-1, 0)} and the 
coefficients are given by 

,(2) _ (2ni + ?i2 + 71 +72 +73 + 2) (2ni + n2 + 71 + 72 + 73 + 3) 



D 



D 



D 



D 



~ (2711 + 2712 + 71 + 72 + 73 + 2)(2?ii + 2712 +71 +72 + 73 + 3) 
(»i + 71 + 72 + l)('^i +72 + 1) 
(2711 + 71 + 72 + l)(2ni +71+72 + 2) 

(2) (»2 + 73 + l)(2?ll +712 + 71 +72 + 73 + 2) 

(O'l) " (2711 + 2712 + 71 + 72 + 73 + 2)(2?ii + 2n2 + 71 + 72 + 73 + 3) 
(2711(711 + 1 + 71 + 72) + (71 + l)(7i + 72)) 
^ (2ni +71 +72 + 2)(2ni +71 +72) 

(2) (»-2 + 73 + 1) (»2 + 73 + 2) 

(-1^2) - (2711 + 2712 + 71 + 72 + 73 + 2) (2ni + 2n2 + 71 +72 + 73 + 3) 

^ 7ii(ni + 71) 

(2?ii + 71 + 72 + l)(2?ii + 71 + 72) 

(2) _ 2(2711 + ^2 + 7i + 72 + 73 + 2)712 

~ (2711 + 2712 + 71 + 72 + 73 + 3)(2ni + 2n2 + 71 + 72 + 73 + 1) 

(7I1 + 1 + 71 + 72)(71l +72 + 1) 

^ (2711 + 71 + 72 + l)(27ii +71 + 72 + 2) 
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^(2) 2(^2 + 73 + l)(2wi + 712 +71+72 + 1) 

(-1^1) ~ {2ni + 2?i2 + 71 + 72 + 73 + 3)(2ni + 2n2 + 7i + 72 + 73 + 1) 

^ ni(ni + 71 ) 

(2^1 + 71 + 72 + l)(2ni + 71 + 72) 

^(2) (2ni + 712 + 71 + 72 + 1)»2 

(o.-i) - (2ni + 2n2 + 71 +72 +73 + 2) (2ni + 2n2 + 71 + 72 + 73 + 1) 
(2^1 (ni + 1 + 71 + 72) + (71 + l)(7i + 72)) 



D 



(2ni + 71 + 72 + 2)(2ni + 71 + 72) 

(2) _ n2(Tl2 - 1) 



D 



(2ni + 2n2 + 71 +72 +73 + 2) (2ni + 2n2 + 71 + 72 + 73 + 1) 
{ni + 7i + 72 + l)(?ti +72 + 1) 
(2711 + 71 + 72 + l)(2ni +71+72 + 2) 

(2) _ (2ni + ^2 + 71 + 72 + l)(2ni + ^2 + 71 + 72) 

~ (2711 + 2^2 + 71 + 72 + 73 + 2)(2ni + 2^2 + 71 + 72 + 73 + 1) 

nijni + 71) 

^ (2771 + 71 + 72 + l)(277i + 71 + 72) ■ 

The operator £" can be written as 

= D^il^iE^, - 1) + i?g;_i)(i?„-i - 1), 

where 

^(1) (2?ii + 712 + 71 + 72 + 73 + 2)(?l2 +73 + 1) 

(0^1) ~ (2711 + 2772 + 71 + 72 + 73 + 2)(2?ii + 2n2 + 71 + 72 + 73 + 3) 

and 

^(1) n2(2ni +712+71+72 + 1) 

(0-1) - ^2ni + 2712 + 71 + 72 + 73 + 2)(27ii + 2712 + 71 + 72 + 73 + 1) ' 

A. 3. Krawtchouk polynomials. When p = 2 equation (|5.19p gives 

K2{n;x;p;N) = - — — ^ kn^{xi;pi; N - 112)^,12 (2^2; ;N - xi 

[~J-^)ni+n2 \ 1-pl 

and the following formulas hold for the generators of Ax 

£^ = pi2;2A:^iV^2 +p2XiVx^Ax2 +pi{xi -N)Ax, + (1 -pi)a:iV^, 
+ p2(a;2 - N)Ax, + (1 - p2)a;2V:,, 

£f = -^(Xl +X2^ N)Ax, + ^-^i— ^X2V,,. 

i — pi i — pi 
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For An we use (|5.25p and we obtain 

i — pi — Pi 

, P2 , 1 -Pl -P2 „ 

+ 1 —(^12 - A^)A„, + — n2V„, 

1 — Pl i — Pl 

Pl(l -pl -p2) ^ ^ plp2\ 

H ("1 - N)A,ai + 1 - Pl + ni^m 

i - pl \ i - pl/ 

£^ = pi(ni +n2- N)An, + (1 - pi)ni V„,. 
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